On The Low Frequency Quasi Periodic Oscillations of X-ray 

Sources 

CM. Zhang 

National Astronomical Observatories, Chinese Academy of Sciences, Beijing-100012, China 

zhangcm@bao .ac.cn 



ASTR1245— 2nd revised— July 2004 
ASTR1245— 3rd revised— Sept 2004 



Received 



accepted 



- 2 - 
ABSTRACT 

Based on the interpretation of the twin kilohertz Quasi Periodic Oscillations 
(kHz QPOs) of X-ray spectra of Low Mass X-Ray Binaries (LMXBs) ascribed 
to the Keplerian and the periastron precession frequencies at the inner disk 
respectively, we ascribe the low frequency(0.1-10 Hz) Quasi Periodic Oscillations 
(LFQPO) and HBO (15-60 Hz QPO for Z sources or Atoll sources) to the 
periastron precession at some outer disk radius. It is assumed that both radii 
are correlated by a scaling factor of 0.4. The conclusions obtained include: 
All QPO frequencies increase with increasing accretion rate. The theoretical 
relations between HBO (LFQPO) frequency and the kHz QPO frequencies are 
similar to the measured empirical formula. 

Subject headings: X-rays: accretion disks — stars: neutron — X-rays: stars 



-3- 



1. Introduction 

The QPO mechanisms of low mass X-ray binary (LMXB) have been paid much 
attention since early 1996 (van der Klis 2000) following the discovery of kHz QPOs by 
the Rossi X-ray Timing Explorer (RXTE). The Z sources (Atoll sources), which are high 
(less) luminous neutron star LMXB (Hasinger & van der Klis 1989), typically show four 
distinct types of QPOs (van der Klis 2000). These are the normal branch oscillation 
(NBO) ~ 5 — 20 Hz, the horizontal branch oscillation (HBO) Uubo — 15 — 60 Hz, and 
the kHz QPOs U2{ui) ~ 200 — 1200 Hz that typically occur in pairs with upper frequency 
1/2 and lower frequency ui. In several Atoll sources, nearly coherent burst oscillations 
Ub-arst — 300 — 600 Hz have also been detected during thermonuclear Type I X-ray bursts. 
These are thought to be at the spin frequency of the neutron star (NS) or double it (see 
e.g., Strohmayer et al. 1996, Strohmayer & Bildsten 2000, Muno 2004). Moreover, the 
low-frequency (0.1 - 10 Hz) QPOs that have been detected in the accreting X-ray black 
hole candidate (BHC), and the QPO properties of BHCs show many similarities with those 
of NS sources and the QPO frequencies of NS and BHC follow a tight and systematic 
correlation over three orders of magnitude in frequency (Psaltis et al. 1999, Belloni et al. 
2002). Nevertheless, the frequencies z/2 and ui, as well as 1/2 and i^hbo; follows very similar 
relations in five Z sources (Psaltis et al. 1998). 

However the frequency separation between the upper and the lower kHz QPOs 
Av = 1/2 — 1/1 is not a constant lower or higher than the inferred NS spin frequency, but in 
some cases decreases systematically with instantaneous M, e.g. Sco X-1 (van der Khs et al. 
1997), 4U1608-52 (Mendez et al. 1998a,b), 4U1735-44 (Ford et al. 1998) and 4U1728-34 
(Mendez and van der Klis 1999). In the latter, the observed coherent burst frequency 364 
Hz is higher than its maximum Ai/ ~ 355 Hz (Mendez and van der Klis 1999). 

Various theoretical models have been proposed to account for the QPO phenomena in 
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X-ray binaries (for a review see e.g., Psaltis 2000). In the early detections of the RXTE, 
the upper kHz QPO z/2 was simply considered to originate from the Keplerian orbital 
frequency at the preferred radius close to the compact object that experiences the inner 
accretion flows, and the lower kHz QPO ui was attributed to the beating of this frequency 
with the stellar spin frequency Ug (Miller et al. 1998). However recently, general relativistic 
effects have been invoked by Stella and Vietri (1999, hereafter SV99) to account for kHz 
QPOs, and these can satisfactorily explain the varing kHz QPOs separation Au. In this 
model of SV99, the twin kHz QPOs are ascribed to the Keplerian frequency uk and 
the periastron frequency of material orbiting the central mass at some disk radius, i.e., 
1^2 = = (M/47r^r^)^/^ and ui = ukW — — 6M/r)^/^] , where r is the Schwarzschild 
coordinate distance and M is the gravitational mass of the source (We set the speed of 
light and the gravitational constant c — G — 1 throughout this paper). It is commonly 
believed that the kHz QPOs are produced close to the innermost stable circular orbit 
(ISCO) or the surface of NS, which provides a probe to detect the accretion flow around 
the non- Newtonian strong gravity region (van der Klis 2000, SV99). 

A theory of epicyclic parametric resonance in a relativistic accretion disk was proposed 
(Abramowicz & Kluzniak 2001, Abramowicz et al. 2003), in which the twin kHz QPOs 
occur at the frequency of meridional oscillation and the radial epicyclic frequency in the 
same orbit. This can explain the frequency ratio 3:2 detected in black hole candidates 
(e.g., XTE J1550-564, Remillard et al. 2002). Although many other viable new ideas 
have also been proposed (e.g.. Wagoner 1999, Psaltis & Norman 2000, Osherovich & 
Titarchuk 1999, Titarchuk et al. 1998, Klein et al. 1996), there has not yet been any model 
satisfactorily explains all observed QPO phenomena. The QPO mechanisms of accreting 
X-ray binaries are still debated and there remain open problems. In the 1980's, HBO 
frequency (z^hbo — 15 — 60Hz) was interpreted as the beat frequency between the Keplerian 
frequency of the magnetosphere-disk and the stellar spin frequency using the standard beat 
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frequency model (BFM) (Alpar and Shaham 1985, Lamb et al. 1985). Later, this was 
considered to be the nodal precession due to the Lense-Thirring effect in the disk (Vietri 
and Stella 1998). Further, in the model of Titarchuk and Wood (2002), HBO (as a low 
frequency lyiow) was ascribed to the magnetoacoustic oscillation frequency, and then a linear 
correlation between this low frequency and the high frequency i^high (identified to be the 
Keplerian frequency) is inferred. 

In this article, we concentrate on the explanation of LFQPO (viow) or HBO of luminous 
NS X-ray sources (for Atoll sources, 15 - 60 Hz QPO is supposed to be due to the same 
mechanism as HBO of Z sources, see e.g., Psaltis et al. 1999), and their relationship to 
the twin kHz QPOs. However we neglect the details of the physical mechanism for QPO 
production at the present time. A summary and conclusions are given in the final section. 

2. The Model and its Results 

The accretion flow of the matter around NS and BH is complicated, especially close to 
the innermost stable circular orbit (ISCO), the radius Risco — 3i?s, where Rs — 2M ~ 3m 
(km) (one solar mass corresponds to about 3 km of the Schwarzschild radius) is the 
Schwarzschild radius calculated using the gravitational mass M, where m = M/Mq is the 
mass in units of solar mass. In reality, the motion of the accreted matter will be influenced 
not only by the gravitational field but also by the magnetic field. For simplicity, we 
assume that the disk matter that exhibits QPO is mainly dominated by the Schwarzschild 
gravitational field with a slightly eccentric orbit (e ~ 0). 

Here we assume that HBO frequency vebo, as well as LFQPO t^qpo, is the periastron 
precession frequency of the accreted orbiting material at some outer disk radius, and the 
twin kHz QPOs are ascribed to the Keplerian frequency and the periastron frequency of 
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the orbiting material at some inner disk radius (as described by SV99). It is also assumed 
that there exists a scaling factor to connect the two radii, which will be determined later 
by the comparison of the QPO data with the model. Therefore, these QPO frequencies are 
conveniently arranged as follows by defining the parameter y = which is the ratio of the 
Schwarzschild radius to the instantaneous inner disk radius. 



u,(y) = (1 - yjl - 3y) X u^iy) , (1) 
i^2{y) ^ i^oV'/' , (2) 



11300 ,^ , 

= {Hz) , (4) 

m 

where is the scaling parameter connecting the radii of the inner disk and the outer disk, 
which we set to be 0.3 ~ 0.4 to give the best consistency when comparing the model with 
the observational data. 

From Eq.(l) and Eq.(2), the parameter y can be expressed as, 

y = -(-)[2-^], y = 0.2{^fK (5) 

^ 3^iy2 V2 1000^ ^ ' 

Considering a maximum value y=l/3 at ISCO from Eq.(l), the parameter y can be 
approximately acquired using a Taylor expansion of (py to the first order in Eq. (3) since the 
maximum value of 0y ~ 0.4/3 ~ 0.13, hence 



(6) 
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At the innermost stable circular orbit (ISCO), the parameter y — ^ and the maximum twin 
kHz QPO frequencies are given by I'Max = ^^1(1/3) = ^2(1/3) = '^^{Hz), corresponding 
to 1200 (Hz) for 1.9 solar masses. The mass of the source can also be expressed by the 
maximum kHz QPO frequency at ISCO, 

2175 

m — . (7) 

^Max 

The relation i^hbo vs. V2 is plotted in FIG.l, together with the measured five Z-source 
samples, and it is seen that the agreement between the model and the observed QPO data 
is quite good for the selected values of the NS mass of about m=2.0 (M©) and oi (j) — 0.4, 
which are the two free parameters in the Eqs.(l-3). 

In FIG. 2, we take the low-frequency QPO vqpo of Atoll-sources, as well as several other 
NS systems and a number of BH binaries for a similar mechanism to HBO of the Z-sources 
(Psaltis et al. 1999), and plot the theoretical curve with respect to lower kHz QPO v^. 
We find remarkably consistent correlations between the theoretical expectations and the 
detected data which extends over nearly three orders of magnitude in QPO frequency. 
However one interesting fact is that the theoretical curves are only weakly related to the 
stellar masses because Vi and Vqpo are all inversely related to the mass parameter in Eqs.(l) 
and (3). Their ratio vi/vqpo — 0~^/^(l — -\/l — 2>y) /(I — ~ 307/) has no direct correlation 
with stellar mass. Further we have for = 0.4, 1 < LOG{ui) — LOG{vqpo) < 1-25, which 
means that the theoretical curves show little variation in the i^qpo vs. vi diagram as a log 
scale even if the mass changes by one order of magnitude. We plot theoretical curves with 
parameter conditions m — S.OM© and = 0.4 (short dashed hne), as well as m = 10. OM© 
and = 0.4 (long dashed line), and find that the theoretical curves are insensitive to the 
mass parameter. However, the mass parameter will influence on the maximum frequency of 
the lower kHz QPO at ISCO because it is inversely related to the mass from Eq.(7). For 
example, if 300 Hz QPO in GR 1655-40 were the maximum frequency, then the mass of this 
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source would be 7.0 solar masses, inferred from Eq.(7). It is also seen that i^hbo increases 
slowly with i/i when ui is higher, as anticipated by Psaltis et al. 1999. 

The weak correlation with source mass in the i^qpo vs. ui diagram seems to hint that 
the QPO phenomena of Atoll sources, Z-sources, as well as other accreting X-ray binaries, 
are intimately related to the specific radii of the disk and the ratio between them. This may 
reflect a common property of the accretion flow around the gravitational sources. However 
the mechanism to account for this is still unclear. 

Prom Eqs.(l), (2) and (3), we can derive the theoretical relations between QPO 
frequencies which are the following, 

^HBO ^ ^Hz) (^)[1 - 0.19(1 - m^r^'] , (8) 



.HBO - 297.30^/^ (Hz) m'^'i:^f'[l + O-150(^)^/^] , (9) 



ly, ~ 297.3 (Hz) m'/'{-^)'/'[l + 0.15(^)^/3] . (10) 

^ ^ 4ooo^ ^ 4ooo^ ^ ^ ^ 

If we set the orbit scaling factor = 0.4, we have the following specific relations. 



i/HBO - 50.6 (Hz) (— )[1 - 0.11(^)2/5] , (11) 
HBO V ^ v^QQn ^500^ ^' ^ ^ 



.HBO - 30.1 (Hz) m^^mi + 0-06(^)^/^] ■ (12) 

As a comparison, we list the empirical relations between HBO and twin kHz QPO 
frequencies (Psaltis et al. 1998, Psaltis et al. 1999), 
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z/HBO ^ (42 ± 3 Hz){ui/500 if^)0-95±o-i6 , (13) 

with 02 ~ 4.6, and 62 ~ 1-8. Therefore a consistency between Eq.(ll) and Eq.(13), as well 
as between Eq.(12) and Eq.(14), is apparently found. Furthermore, we can also obtain the 
NS mass formula represented by the twin kHz QPO frequencies from Eqs.(l-2), 

^500^ HOOO^ ^ 2u2 
According to Eq.(15), we can apply the detected twin kHz QPO data (van der Klis 2000) to 
calculate the stellar mass and obtain the mass average value 1.86 Mq for the six Z-sources 
and eleven Atoll sources, which is close to the measured NS mass 1.78 ± 0.23(Mq) of X-ray 
binary Cygnus X-2 (Orosz & Kuulkers 1999). 

3. Summary and Conclusions 

We have derived a theoretical relation between HBO and ui, as well as HBO and 1^2, 
and found that the detected low-high frequency correlation (Psaltis et al. 1999, Belloni 
et al. 2002) can be interpreted if i/igyj and Vhigh a-re identified as HBO frequency and ui 
respectively. Furthermore, from Eq.(ll), if we set the parameters m=1.8 Mq and Ui = 700 
Hz, then we have the approximate relation ubbo — O.OSui, which is very close to the 
detected empirical relation given in Eq.(13). The correlation between i^hbo a-nd i/i depends 
weakly on the properties of sources, such as mass, magnetic field and the hard surface 
of compact objects, which has been anticipated by the observations (Belloni et al. 2002, 
Titarchuk & Wood 2002). 
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in conclusion it is remarked that the model described here is a simple and rough one. 
Many physical details have been neglected. These include NS spin induced gravitomagnetic 
effect, NS quadruple induced nodal precession (Vietri and Stella 1998), the self-gravity 
of the disk, magnetosphere structure and magnetic axis inclination, the spiral-in effect of 
the accreted matter, the origination and influence of the non-zero eccentricity, etc.. In 
particular, the non-zero eccentricity should have some effects on the QPOs but how this 
would originate is still unknown. We can at least speculate that the motion of the accreted 
matter in the disks might not be described exactly by a free test particle in a circular 
orbit of a purely gravitational field. Consideration of these factors will inform our future 
exploration and understanding of QPOs. Nevertheless, our proposed QPO frequencies are 
all inversely related to the radii of the accreted disks, and the instantaneous accretion disk 
radius is inversely related to the mass accretion rate (see e.g., Shapiro & Teukolsky 1983), 
which suggests that the QPO frequency would be proportional to the mass accretion rate. 
This conclusion is consistent with the QPO detections of Z-sources and Atoll sources as 
well as BHCs (van der Khs 2000). 

Thanks are due to T. Belloni, M. Mendez and D. Psaltis for providing the QPO 
data, and discussions with T.P. Li and W. Zhang are highly appreciated. The author is 
thankful to the anonymous referee for critical and valuable comments, and suggestions for 
the revision of the paper. 
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Fig. 1. — HBO frequency versus the upper kHz QPO frequency for five Z sources of LMXBs 
(Psaltis et al. 1998, 1999 and references therein). Error bars are not plotted for the sake 
of clarity. The model presents a good consistency for the nearly circular orbit of NS mass 
about 2.0 solar mass with the scahng parameter = 0.4. 
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Fig. 2. — The lower kHz QPO frequency versus the low QPO frequency for the Z-sources, 
Atoll sources and other sources (Psaltis et al. 1999). Error bars are not plotted for the 
sake of clarity. The theoretical curves are only weakly related to the NS mass parameter. 
The sohd(dot) line represents the parameter condition m=2.0 and (f) — 0.4 (0 = 0.33); the 
short(long) dashed hne represents the parameter condition m=5. 0(10.0) Mq and = 0.4. 
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